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NOTATION 


Cc Chord of control plane section 
Fe 
c. = Lift coefficient of control plane section 
puéc 
2 
D Diameter of the submerged body 
Ey Exponential Integral 
Fx Axial Force 
Ey Vertical force 
: Px 
Fy? Non-dimensionalized axial force 
x 212 
pusL 
2 
| P, 
os Ps Non-dimensionalized verticai force 
2,2 
s puch 
a 2 
2 
=| u 
4 ty Vai Froud# number of the sulmersible 
a 
4 U 
ie le Tec broude number of the controi plane section 
g Gravitational accelerat.on 
G fhree-dimensional Green function 
G, Two-dimensional Green function 
h Depth to the axis of revolution of the aubmerged body 
h, Depth of submergence to control plane section at quarter chord 
3 ho Depth of submergence to control plane section at mid chord 
: i = ¢(-i) Imaginary Unit 
| L Length of the submersible 
Hy Moment about the oy~sxis 
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q Unit vector of the body surface drawn into the fluid 
rs) Freasure 
r ‘ Position vector of the center of gravity 
U Speed 
p Bh Water denaity 
a angle of attack of the control plane section 
i] Pitch angle cf the bare huil 
r Vortex strength 
Go Source or sink strength 
6 Three-dimenstonel velocity potential 
& Two-dimensional velocity potential 
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The forcea and moments acting on a 
submersible are computed when it is 
beneath the free surface or near a wall. 
The method used in this report is based 
on potential theory. In the computation 
of forces and moments on the hull, a 
three-dimensionsal sethod is applied. The 
free surface condition tg linearized and 
the body boundary condition 4a exact. 
The body surface is discretized with 
surface elements and the singularity of 
source and sink is distributed on thea. 
The strengths of the unkown sources and 
sinks are determined through the bedy 
boundary condition. A two-dimensional 
method ts applied £9 compute forces and 
moments of theyceatrot surfaces. The 
boundary conditions arel's'ame ee those for 
three-dimensional case. Computed forces 
and momente of the control planes are 
added to those of the bere hull. To 
include the interference effect of the 
hull on the control planes, the flow 
velocity at each control plane is 
computed with existence of the hull. 
There are some discrepancies between 
computed results and experimental data 
because of the effect of viscosity. The 
overali trends in the computed results 
are same as those of the experiments. —— 
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INTRODUCTION 


There is a need to develop an analytical method to compute the forces and 
moments acting on a submersible when it ite beneath the free surface or near a 
wall. The methodologies presently available to compute these hydrodynamic 


forces and moments are limited tn epplication. Pond! # computed the momenta 


developed on a Rankine ovoid by using the method of axial distribution of 


sources and sinks. McCreight. improved the method used by Pond by distributing 
the dipoles and computed the vertical force. Their results are in good agreement 
with the experimental data. However, the methods used by Pond and McCreight are 
valid only for zero pitch angle. 


The present method computes the vertical and longitudinal forces and nitch 
moment acting on a submersibie with and without control planes when the 
submersible moves beneath the free surface or near a wall. This method is more 
accurate mathematically than that of Pond: « Furthermore, the effect of pitch 
angle ig incorporated in the computation. For the computation of ferces and 
moment acting on the bare hull, the so-called panel method is applied. The body 
surface is discretized with many quadrilateral planes, and sources and sinks are 
distributed on these surface elements. This method can be applied to a 
submerged body of arbitrary shapee To compute the forces and moment of the 
control planes, the method developed by Glesing and Smith’ is used. 
Two-dimensional sources and sinks are distributed around the sections of a 
control plane, end a vortex located in the middle of each section is introduced 
to compute lift force. The inlet velocity to the control planes is computed at 
the tips of the sections of the control plane to include the interaction effort 


between the bare huli and control planes. 


*References are given on page 45. 
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VELOCITY POTENTIAL 
HULL 


The coordinate system, oxyz, moves at 2 speed U, which fe the mean forward 
speed of the submersible along the positive ox-axis (see Figure 1). The posi- 
tive oz~axis {6 always directed verticelly upwards. The origin O is located 
above the center of gravity of the submersible. The oxy-plane {6 the plane of 


the undisturbed free eurface. 


We asaume initially that the submersible and coordinate system are sta-~ 
tionary and that the fluid around the submersible moves toward the negative ox- 
axis wit uniform epecd U. Then, the total velocity potential for the bare hull 


cen be expressed. by 
O(x,y,z) = -Ux + (x,y,z) (1) 


where 6 is disturbance velocity potential due to the submersible. The distur- 


bance velocity putential satisfies the following conditions: 


1. Leplace equation in the fluid domain 


2 2 2 
Be oe, 28g (2) 
6x oy oz 


7h The linearized free-surface condition 


2 
: fe + es = 0 (3) 
- dx (Oz 


3. The body boundary condition 


Ot yh (4) 
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4, The radiation condition: the disturbance vanishes sufficiently fast 
far shead of the ship. 


cy 2 2 
fia RCO eet Oe) 10 (5) 
x72 


Re oe n y2) 2 


5e The bottoa condition: the normal velocity at the bottom is zero. 
\ 
b2(%,y02) =0 asz+-e (6) 


In equation (3), the censtant ie given as 


U : 


and in Equation (4), the unit vector which is orfented normal to the fluid is 
given by 


(ny> Ros n,) =n 


4 
| 
: 
- 
; 
: 
# 
a] 
< 
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— (8) 
(ny> Ne» no) =a Fe za 
»> 
where ive is the position vector of the center of gravity. 


4 
The solution of Equation (2) is given by Brard in an integral form as 
1 : 
(x,y,z) = --~ jf G(P,Q)0{Q)48(Q) (9) 
an S 
where P(x,y,2) is the field point, Q(xgr¥or2g) the source point, S the wetted 
surface of the body, o the unknown strength of sources and sinks distributed on 


the body surface, and G is the Green function which ie given by Wehausen end 
é A 
Laitone” as 
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w 2 
G(P,Q) = Ae Ko sec 6 dO f 2 du (10) 
y a es CRP Unik Be CuO 
where 
: 2 2 2 2 
Ce (Xe KGa try yg) alt GZ) eazy) (11) 
5 2 2 2 2 
Bie CEE Gs tC yee ty gy) hil (zu tanz Gi) (12) 
. and 
@ = (x = x9) cos® + (y - yo) sind . (13) 


Only the real part of Equation (10) is used in later computation. 


The strength of the sources and sinks, ¢o, in Equation (9) can be found hy 


substitution of Equation ($) into Equation (4) to obtain 


zal ond 
a 


aa {f oCP.QdoCayds(a) * Un, | (a4) 


The solution of Equation (14) is only feasible with the help of a numerical 


3 procedure which will be given later. 


The force acting on the body is expressed as 


i F=-ffpa ds (15) 
S 


where p is the pressure around the body and is linearized from Bernoulli's 


equation as 


ies 
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p=p U 4, (16) 
where p fe the density of the flutd. 
CONTROL PLANE 


The velocity potential for the control plane, whose coordinate system is 


shown in Figure 2, is expressed in the two-dimersional domain as 
d(x,z) = -Ux + 4(x,2) (17) 


where ¢@ is the disturbance velocity potential due to a section of the controi 


plane. This disturbance potential satisfies the following conditions. 


le Leplace equation in the fluid domain 


ra aa 
oe, Seo (18) 
6x éz 


2e The linearized free-surface condition 
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oe tks be o (19) 
Ox oz 


3. The body boundary condition 


én I 


4. Kutta condition: the velocities at the tralling edge elements, one at 


the top cf the surface ana one at the bottom, are equal. 


In addition, the disturbance potential should satisfy radiation and bottom 
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conditions similar tc Equations (5) and (6). 


eee cere eee =e ea a a om 


Oh na 
LOG eye 
\ iti 


hob adhneh pas iunmnad, 


Bree tt ts 


bal ya eset 


Lec ACRE 


Rte eae Pe on 


Beg anche ei SNe 


3 
The solution of Equation (13) is given by Giesing and Smith . The velocity 
potential function is divided into three parts as 


o(x,2) = o) + Mo, + %,)- (21) 


The velocity potential v) is due to sources and sinks which are distributed on 


each section of the control plane. Vo is the velocity potential due to a vortex 


of unit strength which is located at the inside of a section. The vortex is 


introduced to calculate the lift force around the gection. Ys the last term of 


Equation (21) is the velocity potential due to sources and sinks of unit 


strength distributed on the section. This velocity potential cancels the normal 


velocity generated by $5. If we substitute Equation (21) into Equation (20), 
the following conditions are given 


% 
le vn, (22) 
on ; ( 
and 
Ob2 043 
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These three velocity potentials are expressed as 
b.(P) = 4 f 6.(q)6.(p,q)d1(q) for i=] and 3 (24) 
i 27 ac? 2 


and 
oP) = Re{ i. in{ x-a + i(z=b) | = son [ (x-a) + i(ztb) | 
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where o4 fe the two-dimengional strength of sources and sinks, (a,b) the 
location of the vortex, q(x9,%9) & two-dimensional source point, and Go is the 


two-dimensional Green function given by Wehausen and Laitone> as 


Go(p.q) = Re {In [x-xg + i(z-z9)] + In[x-x, + 1(2+z9) | 


(26) 
a e7ik[ x-xo + 1(z+z,) | dk Bo Eriko x=XG + (2429) ], 
+ 2pv J k~ ky 


In Equations (25) and (26), Re denotes the real part of a complex quant/ty and 
py the principal-value integral. The vortex strength I in Equation (21) will be 
determined with the Kutta condition. The pressure around each section of the 
control plane can be obtained from Equation (16) by substituting 4, for o° The 
force and moment acting on the section can be computed with an equation of the 


same form as Equation (15). 


By substituting Equation (24) into either Equation (22) or Equation (23), 
the unknown etrength of the sources and sinks for the two-dimensional case can 


be determined from the following equation 


ne J 04(4)S,(p,q)d1(q) = Un for i= } : (27) 
Cc 


26 On 1 


eYa 
we —--—~ for {=3, 
on 


When the submersible is moving near a wall, the linear free-surface 
conditions, Equations (3) and (19) are no longer valid. Instead, we should use 


the wall boundary conditions given by 
oo 0) at z= 0 (28) 
Co) 


and 


AES 


exeneab 


fe Ae eon (29) 
dz J 


All other boundary conditions are the eame ag those described in Equations (2), 
(4) - (6), (18) and (20). The solutions in this cage have the same form as 
Equations (9) and (21), with different Green functions. The Green functions can 


be expressed by application of the method of images (see Reference 5) as 
\ 


G(P,Q) -+ + = | (30) 


for the case of bare hull and 
Go(P.q) = Re(In[x-xq + 1(z-z5)]} + In [x-xg + i(ztzg)]} (31) 


for the case of control plane. The velocity potential due to a unit vortex 


located at (a,b) should be 


be) = Re (— Inf (x-a) + i(z-b)] - = In [(x-a) + 1(z+b) ]} (32) 


To find the velocity potentials for the wall condition, Equations (19), (25) and 
(26) should be replaced by Equations (30), (32) and (31), respectively. 
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NUMERICAL PROCEDURE 


To determine the unknown strength of the sources and sinks in Equation (14), 
we first discretize the wetted surface of the body S with many quadrilateral 
elements. Figure 3 shows the discretized surfaces of a spheroid. Furthermore, 
we assume that the unknown strength GE the sources and sinks, o(Q), is constant 


on each surface element. ‘Then Equation (14) can be written as follows 


N 


a) 2 
ie sh. 04(Q4) Sq CCPL 994) AS, a “Un, for £ 1, 2, coos N (33) 


i 


Once the derivations of the Green function, G, are evaluated sumerically, 


Equation (33) can easily be solved for Oye The numerical evaluation of the 


E Green function and its derivatives is the most difficult part, and consumes a 


lot of computer time when running the program. 


The normal derivative of the Green function in Equation (33) can be 


expressed as 


eG 2G G G 
ae 80 26 oe , G4) 
on 6x Oye he Oz 


The task is to find three derivatives of the Green function. To do this, we 


first derive the Grean function suitable to numerical evaluation. The Green 


function, Equation (10), can be rewritten as 


Tt 
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where 


O, = (x-X9) cos @ + (y-y,) sin 0 2G 
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(x-%9) cos @ - (y~y,) sin @ 
(36) 
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The notetion under two of the integrals in Equation (35), namely, Ly or Ly» 
represents the integral paths shown in Figure 4. To derive the derivatives of 
the first two terms of Equation (35), the method which Hess and Smith’ developed 
for unbounded flow can be used. Kong and Paulling’ heave applied this method for 
the computation of motions ef a body in wavea, and their paper can be examined 
for the details. Since Reference & can be examined for these detaile, we will 


present here the derivation of the last two terms of Equation (35). 


We write Equation (35) as follows 


i 1 
wad e 3 
2 @ .72U,_iu,u WU 
Ge-t f uae { f2—& te) an 
o o Sees 
(L)) (38) 

é Zu ~iwyu He gi®aly 
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We now change the varfable for the integral path, Li» as follows: 


v=(u-ug) (2-10) with & © %4 or O.; 
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then, the integral path becomes that shown in Figure 5. The integral path, Ly» 


ic now traneferred onto the path C, in Figure 5 aes 


© Zu (0s — _ 4% -v 
{= CET Cy me a etl | eS (39) 
(L)) 


The integral over c can be easily evaluated with the method of residue as 
e~ e~ Ini; <0 
J <—av= son dv. +. | ha re (40) 


with the introduction of the exponential integral 


|= dv = BE, [-uple = iw)], (41) 
v 
“Uo(z = 10) 


the integrai over Cc, becomes 


-v 4 one 
f= dy elas - ioe t |. 2 
Cc 
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By substituting Equation (37) into Equation (33), the integral over L, ie given 
by 


“ZU, , iW.u Lou = ~~ 
Qn (Omer rete erica) ~ evTtol® = iw.) _ ~ _ 2ni 
f u = Us rae are {E, [-19(z 104) } 0 } 
(Ly) 
(43) . 


tol? - iw) {z,[-uolz - ia.) + (ease Oy, W. < 0 


0 We, We > O 


For the integral path, Lys we change the variable as follows: 


vefuc ig) (z + iu) with we O4 or On3 


then, the new integral path is given in Figure 6 With this new variable, the - 
integral path, Los is now Cie Enploying the sane procedure ae used previously 


to evaluate the integral over Ly, the integral over Ly becomes 


oe -fu 104u ~i0u = no 

a Ce e ) un(z + iw) - ~ and 
BN nn hy ° ef ee + 

f Sane du= 6 {z [up (z iv,) | 0 } 
e u 

(Ly) 


-u,(z + iw oN aa Int) Wy, On 
els Ug t te) te, [uC Aye aa Wy, W. <0 (44) 


0 Oye w. > 0 
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Substitution of Equation (43) and (44) into Equation (38), allows us finally to 


express G3 in a convenient fore for numerical evaluation ag 


a/2 


el 
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f 
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7 en YolZ 


uo (Zz 
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If we let 


—u,(z - iw 
Ue of +) 


ger noe ~ iw) 


- ~- dw. 
wpe Ug (Zz U3.) 


— ~ 
-u,(z - tu) 
of 


uae de eke o tO) fg [ug (z fh 104) ] 5 any 
0 


= fuw_) {B,[ug(z - 102] % cat 
0 


+ iw.) {E,[-ug(z + Poot zat 
0 


+ 1w_) fe, [-ug(z + tools ate oe te <0 
Wr, W > 0 

BE [-uo(z - tuy)] = R + iI 

x#C+is 

zB [ugl2 = tu.)] = R + 41 


eC +) 14S 
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(45) 
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we can express Gy as 


n/2 


The derivatives of G, can be eaally 


where 


6G 


—e» ft u cosd G, 
x 


3G, ~ 
—e-= iu sind G, sign(u) 
ay 


sign (9) = ] for 


= -1 


Ee 
a 
= 
+ 


Et 
& 

& 

i 


37-2 [r+ CDser u(r) 8, to To (47) 


B, G > 0 


expreesed from Equation (38) 


(48) 


(49) 


(50) 


(51) 
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With substitution of Equation (45) into Equations (48), (49) and (50), the deri- 


vetives of Sy can be expressed in final form as 


® 
éG ‘ ~ 
3 2 ( w = 2 
oar a = a } Ug cosd {- PSR Scc ana a) ¢ ae 


é 2. ~2 
Sis 
Pe eye ar (2%) Ghaa (52) 
=-2 ~2 9 


W 
+ 
€ 


He 
a6, : 


Hane i Wo sind {- “+t - 1 + Gaxe 
dy m of a5 5 re) 
0 came ty CE 
Ue 
+ Coy ty Gee) Chae (53) 
~ 0 
Zu sty 
n/2 
@ 
oG se 
3 ao Zz Zz 25 
bz x J Dodltrmeest Ret (g)5 
> + 
0 meses 
an D) os 
$ ——-—— + B+ (OT) Shae (54) 
oa oy? 0 
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In these equations, tha upper value is taken when &, or & 43 larger than and 
equal to 0 and the lower value is taken when 8B, or OB fie negative. Inglis and 
Price’ have derived the Green function for a translating and puleating aource 


aimilar to G3 with the change of variable and with the change of integral path. 


When the argument, muy (241i) or -u,(z~13) fe small, the exponential integral 


can be eapressed in a series form (sce Reference 10) as 


n 


B(2) ey eine J UES fargz][<x« (55) 
a= is) : 


where z is any coaplex nuaber and y # 0.5772157 is Euler's constant. When the 
argument becomes very large, the exponentiel integral is evaluated by the method 
of toda.2! For a large cowplex nuaber s = x + fy, the exponential integral 


multiplied with exponential function is expressed as 


z 
e EB, (z) = iy - it, (56) 


where 


SU et ut nO tod : (57) 


Is i) Oe ey (58) 
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For the evaluation of q and Ln we use the following approximation 


@e 
I © j e f(t) datz J ve eG) (59) 
0 
fe (a) (n) 
@re Fy are the zeres of the Laguerre polynomial and ry are the 


(5) 


corresponding Christoffel numbers. In thie report, n fe taken as 5, and Ry 
and x? are given in Table I. 


Table I: 
Zeros of Laguerre Polynenials and 
: Christoffel Numbers 


(5) 

Ky yee 
0.2635603 ee 0. 5217556 
1.413603! 1.3986668 
3.5964258 0.0759424 
7.0858100 0.0036i17 

12.6408008 0.0000234 


The exponential functisne tn Equations (47), (52), (53) and (54) are highly 
oscillatory when the resl pert of the complex argument {s smali and the imagi- 
nary part is large. One of complex arguments in Equation (46) can be expressed 


&é 


mu (Z = 10,) = - ko sec” @ + [z + iR cos(@ + 8)] (60) 


13 


ae 


=A iy 
‘ 


dh 


fea ates 


Tk ath gL aS a a a tg a 


(Si es BT i Vt at 


} 
: } 
\ 
3 
4 
4 
4 
; 
3 
i 
‘ 
{ 
| 
4 
A 


Brad ted Dir al a te nr are 


Si eS 


OL SE en le etl ae eS 


with 


=) Youn 0 
6p = tan ‘E> iy 


When e, - 8, is emall, the integral between 6, and 6, can be approximated as 


1 2 
ee 2 
I. ug(C + 18)d0 = ic sec*® « exp{-k, sec 8[z + £R cos(6 + B)]} (61) 
%) Mt 
o2 


= § (A + 4B) exp[i(aé + b) Jaa 


where A, B, a, and b are congetant between 8) and Soe The last expression in 
Equation (61) can be analytically evaluated. Once the normal derivative of the 
Green function, which contains Equation (61), 46 numerically evaluated, Equation 


(33) can be expressed as a system of linear equations and o, can be solved by 


j 


the method of Gauss elimination. 


For the computation of forces and moment of a control plane, we discretize 
each section of the control plane with straight segments as shown in Figure 7. 
The unknown strength ef the two-dimensfonal sources and sinks on each segment ‘s 
assumed to be constant. The subscript 1 for o, in Equation (27) will be dropped 


from row on to avold confusion asseciated with that of segments. o represents 
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each of Sy and 4 in the following equations. Equation (27) can be approximated 


MN 
1 2 : 


Se for £ = 1,2,.00,M 


the normal derivative of the two-dimensional Green function is 


aG2 0s: 8Ga dG 
dn Ox < dee 53 : (63) 


We use the complex varisbles in derivation of derivatives of the Green function 


suitable for nuserical evaluation. If we introduce the following complex 


variables 
C=xtiz and C, 7 Xo + tZo, (64) 
Go(p.q) = Re{F(C,Co)} =. (65) 
where 
Ee . ewik(C = Co) 
F(C,C) = In(C - C,) + In(C - Cy) + pv [ eeu. dk (66) 


0 


Sante Oo %o) 


In Equation (66), C, ie the complex conjugate of Co. Tnen, the two derivatives 


in Equation (63) can be ewxpreased az functions of PF 


O64 


ox 


Rete) ; (67) 
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The derivative of Equation (66) with respect to C is 


erik(C = Gy) (69) 
a. _1_,_1_.t,,[' dh 

8c C= ts C= Co J k = ky 

0 


-i k,(Cc - C,) 
- 2% ky e ‘4 g 


Integrating the last equation with respect to a4 (q4) we obtain 


{ a6 (qq) = -(cos ay ~ 2 sin a;)[1n Ts4] + 1054) - 1n ie 10] 
- (cos o + 4 sin a,)[1n Tis) + feel - in er o 16, | 
kCzi sa cee) 
ee . (cos 8, + i sin a.) {pv — e ae [cos k(x - Xo 
0 


o 


kG + Z9;) 


[cos k(x - X9,) 
k= Pe 


a 
- i sin k(x - %oj4))) 7 py j 
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- i sin k(x - X04) ]} 


ko(z + Foye 


+ 2ni (cos a4 + i sin ay){e coa ko (x - Xo) 


- i sin k(x - Xoi44)] 


Z kote + Zoye1 


cos kyo(x - X04) - i sin k,(x - Xo) J} (70) 


the new variable and notation are given in Figure & The details for derivation 


of Equation (70) are given by References 3 and 12. 


The numerical procedure for determining o, and [ is as follows: First, we 
evaluete the normal derivative of the two-dimensional Green function which 
enables us to solve Equation (62) for Tye Next, we calculate the nornsl 
velocities at segments due to a unit vortex located at the center of the 
gection. To eliminate these normal velocities, we distribute Ba at each segment 
of the section and determine oy with Equation (23). We finally solve for the 
vortex strength [T in Equation (21) using the Kutta condition. 


The computational procedure of forces and moments is ag followe: The forces 
and moments of the bare hull are first computed. This means it is assumed no 
interference of control planes to the bare hulle To include the interference 
effect of the bare hull toa the control planes, the flow velocity, U, in Figure 2 
is different from U in Figure 1. A control plane is cut at three different 
spanwise locations. At the leading edge of each location (or section), U is 
computed. This velocity is different from the steady forward speed, U, in 
Figure 1. ‘Two-dimensional forces and monents are computed at each section; and 
these are summed along the spanwise direction. Finally, these summed forces and 


moments are added to those of the bare hull. 
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“NUMERICAL RESULTS 


In this section we will calculate the forces and moments acting on a subwer- 
sible moving under the free surface or under ice. Numerical results based on 
the present method will be compared with results obtained from earlier theoreti- 
cal methods and experiments. First, forces acting independently on the bare 
hull and control plane will be compared; then, results for the combined hull and 


Plane will be compared. ' 


Figure 9 shows the forces and moments on a Rankine ovoid having a length of 
4 ft (1.22m). The ratio of length to diameter (D) is 10.5. The Rankine ovoid 
is created by distributing sources and sinks along a line which is assumed to be 
in a uniform stream in an unbounded fluid (see Reference 12). The computed ver- 


tical forces (F') are generally somewhat larger than the experimental data for 
Zz 


smaller submergence (h) and smaller for larger submergence. The computed 


twoments (M') are generally smaller than the experimental data for all sub- 
y 


mergences when the Froude number is larger than 0.5. However, the overall trend 
of the computed results is similar to that of the measurements. Results calcu- 
lated by the present method are in good agreement with vaines obtained from pre- 
vious analytical methods. Figure 10 shows the same results plotted as a 

function of submergence. As the submergence becomes larger than three times the 


diameter, the force and moments decrease rapidly for all Froude numbers. 


Figures 11, 12 and 13 show gimilar data for a spherold whose ratio of length 
to diameter is 7. The findings here are somewhat different from those discussed 
previously for the Rankine ovoid (see Figure 9). When pitch angle ts zero (see 
Figure 11) the experimental data for the epheroid extracted from Reference 13 
are larger than the computed values at Froude numbers less than 0.5. The 
agreement between computation and experiment is generally good when the 
submergence is equal to or larger than the diameter: when the ratio of 


submergence to diameter is 0.75, there is substantial difference. When the 
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pitch angle ig 2.5 degrees (bow up), the experimental results are larger than 
those from computation for essentially all Froude numbers (see Figure 12). As 
shown in Figure 13, negative pitch angle of 2.5 degrees (bow down), the 
prediction and measurement of vertical force are, for the most part, in better 
agreement than those for zero or positive pitch angle. The computed moments are . 
larger than the experimental results for all submergences and Froude numbers. 


These discrepanciey for non-zero pitch angle might be caused by the viscosity of 


the fluid and also vortex shedding from the body or the cross flow effect. 


Figure 14 shows the results of lift for a hydrofoil whose section has the 
shape of NACA 4412. The reevlts shown in this figure are those for a two- 
dimensional section. The agreement between computation based on the present 
method and experimental data from Reference 14 is good when the angle of attack 
is smaller than 4 degrees (see bottom figure). As the angle of attack 
incr2ases, the computed results become larger than the measurements particularly 
at the lower Froude number. It is of interest to observe that for this sub- 
mergence the lift coefficient at the higher Froude number is smaller than that 
at the lower Froude number. When the angle of attack is 10 degrees and the 
Froude number ig 0.922 (see top figure) there is substantial discrepancy between 
the computed and experimental results as the submergence increases. The results 
computed by the present method are further compared with the method developed by 
Wadlin and Christopher in Reference 17, Equation (15). ‘Whe calculations were 
made for an aspect ratio of 10 to stimulate a two-dimensional lift coefficient. 
It should be pointed out that Equation (15) of Reference 17 is strictly 
applicable for aspect ratio of 0.125 to 10. Their results are almost half of 
those by the pregent method. Compared with the experiment, the results of 


Wadlin and Christopher show better agreement than those by the present method. 


Figure 15 shows the results of lift for a hydrofcil whose section has the 
shape of NACA 644010. The aspect ratio of this hydrofoil ig 6. As shown ‘n the 
bottom figure, for small angle of attack there are small differences between 

computation and the experimental results taken from Reference 15. The differen- 


ces increase as the angle of attack increases. The reason for these discrepancies 
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may be that firet, the effect of viscosity is not included 4n the computation, 
and second, the test was done with turbulence stimulation. The top figure 
indicatea that for an angle of attack of 4 degrees, the computed lift 
coefficient ig about twice the measured value when the submergence is large. 
When the ratio of submergence to chord is smaller than 0.5, the agreement 
between computation and experiment fig good. For small submergencea, the lift is 
more affected by the free surface than viscosity, and when the submergence is 
large, the reverae is true. On the other hand, compared with the experiment, 
the results of Wadlin and Christopher show better egreement than chose by the 


present method when the submergence ratio is larger than 2. 


When the location of the gail plane or stern plane of a submersible is at a 
depth of two or three times of chord, a better lift computation can be expected 
with the method developed in Reference 17. 


Figure 16 showa forces and moment for Model 4621 with and without Beer 
planes at deep submergencee The computed and experimental axial forces are 
fairly steady at different trim angles. The vertical forces for the bare hull 
are computed to be significantly less than those of experiments However, the 
results of vertical forces with sternplane agree very well with those of experi- 
ment except at a=12°. The reason for the good agreement for the case with 
sternplane igs that the vertical force {in this case lift) of the control plane 
atone is over-estimated at deep submergence as shown in Figures 14 and i5; and 
this over-estimation is compensated with the under-estimation of the vertical 
forces of base hull. The moments are computed to be larger than the experimen- 


tal values. 


Figure 17 shows the results of vertical force of a spheroid with L/D=7 near 
awall. The computed results are compared with the results of Newman. 
Newman's method was developed using slender body theory and with the assumption 
of L/D >> i and h/D << 1. The vertical forces computed by the present method 
are smaller than these computed by Newman. It is unkown which method is more 


accurate. 
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CONCLUSIONS. AND RECOMMENDATIONS 


The forces and moment acting on a submersible moving beneath the free sur- 
face or near a well with a pitch angle for the bare hull or a control surface 
deflection have been computed using potential flow theory. There is good 
agreement between calculation and experiment for the bare hull when the ratio of 
submergence depth to the diameter is larger than 1 and the pitch angle is zero. 
For the control planes alone, the good agreement between calculation and 
experiment is achieved when the ratio of submergence depth to chord is small. 
When this ratio is large, the force due to viscosity is dominant and the present 


method generally computes forces larger than those of experiment. | 
Frem the present study, the following conclusions may be drawn: 


1. The numerical evaluation of the Green function needs to be improved 
when the ratio of submergence depth to diameter is smalier than 1.0. 
Near the free surface, the Green function has oscillatory behavior and 


the exponential function decays slowly. 


26 For nonzero angle of attack, there is a cross flow or vortex shedding. 
The effect of cross flow or vortex shedding should be included in the 


computation with the help of viscous flow analysis. 


3. When the ratio of submergence depth to the mean chord of a control 
plane is large, the lift is dominated by the contribution of viscosity. 
The viscous effect on the lift should be tncotporated in the future 


computation. 


4 To improve the results of potential theory, an empirical or analytical 


approach should te developed with inciusion of viscous effecte 
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Figure 2 - Coordinate System for a Control ‘Plane 
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aUBTg [OI]U0D B Jo WOTIDIg JO VOFILZTIeI19STq - 4 VAN3FY 


queugas 


34 


2 i045 


“yO°O1 


sz 


om 


2881 ‘Nur 41 


LL01d-80° 


‘ 


BSS 


O°5 uldSsid 


jth segment 
1541 541? 26541) 
\ PS . 
a 
95 045 ?705) 
Zh -z 
a, = tan io 21) 
*55+1 oj 
CoH a A) : 2: 
: Bis eS Corey + (2-2, ) 
e 2 LD. 2 
3 Seige edge t Geox aD Cae) 
=) Zaz 
= 6. = ran”! (24) 
_ j xX-xX 
4 Oj: 
=a a Aza 
3 8 Scant seecorrl 
oe aatae ; x-x 
= ott 
4 
wot 
fang 
1 
4 
4 
a: 
S| 
4 
4 Figure 8 ~ Netations for Equation (70) 
4 
af 35 


AMAR 
if \ Y i ni) 
ae 1 Miter a3 
A 1 I pay eS Ati 
ed a He 
| 
; ; ri Nh rey, i i] 
u ea vk i i 
if 5 al (een anon 
j ; ; ; : i! hee 
‘ i ni hl ; of 
i i on Meet 
} i 
¥ i fc 
hae! : ee 7 ; 
iv 
i , 
| j 
t Ha, ste plement See nia, ark 
; ‘ a teen j rant 
i iM y 


ae ey Re y ™ RRR cng ene yale ayn tia pe LP Ln 


f Vy ea 
pile eh 


Pa ea een 3, ce rad 


Mi © primal cgnasiet . 
if 1 4 pik: 1 . 
Vi : ~ : ; f 
‘ is i vt. hay : % i 
Vay aa | Poa 
f ’ y foe 
i Te fo a ; p ys 
i ; f : ran i Vip ry H 
Pt a RY 
saninaipaatli ree a 
i : 1 F 
: : Pavan i 
v 
iy 106 
1 tn) 4 
Q bee ‘ ; ; : y * MNCs 
; { \ ) At SH nN 
fed DeckaNAby seas eam 
: enitey hk oe 
a tty 
1 het 
Saye be 
iy f ; ie mt 
i if ins Aj i isis Renata 
(oul peeme ‘ I ey 
Wty : Hy rare Eto 


Present Test Calc. 
(Ref. 1) (Ref. 2) 


OOO +t 


poo xxx 


Figure 9 - Force and Moment on 4 Rankine Ovoid with L/D = 10.5 


at Different Submergences 
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Figure 10 - Force and Momcnt on a Rankine Ovcid with L/D =.10.5 


at Different Froude Numbers 
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Figure 11 - Force and Moment on a Spheroid with L/D = 7 


when 6 = 0 
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igure 12 - Force and Moment on a Spheroid with L/D = 7 


when 0 = 2.5° (Bow up) 
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Figure 13 - Ferce and Moment on a Spheroid with L/D = 7 


when @ = -2,5° (Bow down) 
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Figure 14 - Lift of a Hydrofoil of NACA 4412 Shape 
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42 


/ 5 
be 
try 
yo: 
Hang 
i 
i 
i } 
‘ 


et yh, 
NDie hike 


TH 


ie 


i 


colina Avutistoinaitelliserpmeranianine ppheyerape, eg 


; ‘ hes] 


vay 
i 
My 


i i" 


A dnp etemahomdrng ig eben oven 


oe t Brae 
P 1 Det ay 
% 
y t 
sii Y 
; A 


2.5 
xl0~ Bare hull | 


Present 5 
Method Test (Ref. 16) 


Q00 
opo00 
+++ 


Bare hull + Stern plane 3A 


Figure 16 - Forces and Moment on Model 4621 at Deep Submergence 
(h/D = 4,89) 


43 


4 2! wai ua 


a aetly a ee Me | 


t 


a etl ake eo IM Tc a Ne ae Sak ae 


ge ry 


Present Method 


ps 


a_—=-—e= Calculation (Ref. 17) 


4.0 


2.0 


0.5 1.0 is 2.0 n/D 


Verney Soars Saher agp aah nrriineh 
cal PE Vit aS Ed a 


1? ‘ 
ely OE 


Vie! RESP doe oe 


os ne SS ete le es a Ne 


teal Force on a Spheroid with L/D = 7 Near a Wall 


aE eis pa at eaters cl 


See a te Na ela ae 


' 1 Nee 
in 
au 
y Jalan enya 
; efi Bare oun 
nan aA: abiding 
i Lt {i fi Nes Paty i ite 
q " hie rte ay aie Moi 
{ 1 " 
i 
. it 
as yay - 
; umes : ‘ Ph 2 Peg aiee nano 
7 vane 4 fy 7 ak ew ae OM , Ts eet re phe 
Fitch nhac aaa RL eke p ti Lickin ei fl onan Merce ae TRCETONR na Rene atayh 
a Puy i Ee pp ie q 
} J ny uh 
ape ny I} 
i sore 
Ls ’ ; 
‘ e y / 
i 
7 ( 
y 
MN crn tM ; 
} SLT Aes By : 
it Scum 8 
; if 
" LGM | kt RR a 


l. 


26 


3. 


4 


56 


7. 


Se 


REFERENCES 


Pond, H.L., “The Moment Acting on a Rankine Ovoid Moving Under a Free 
Surface,” Journal of Ship Research, March 1959. 


McCreight, WeR., “Heave Force and Pitching Moment on a Submerged Body of 
Revolution in Finite Depth,” Massachusetts Institute of Technology, 


- Department of Naval Architecture and Marine Engineering, Keport No. 70-7, 


June 1970. 


Giesing, J. P. and A.H.O. Smith, “Potential Flow about Two-dimensional 
Hydrofoils,” Journal of Fluid Mechanics, vol. 28, pp. 113-129, 1967. 


Brard, R., “The Representation of a Given Ship Form by Singularity 
Distributions When the Boundary Condition on the Free Surface is 
Linearized,” Journal of Ship Research, Vol. 16, No. 1, March 1972. 


Wehausen, J.V. and E.V. Laitone, “Surface Waves,” in Encyclopedia of 


Physics, vol. IX, Springer-Verlag, Berlin, 1960. 


ov 


Milne-Thomson, LeM., “Theoretical Hydrodynamics,” The McMillan Company, 


New York, 1960. ‘ 
Heas, JeLe and AeM.O-e Smith, “Calculation of Non-Lifting Potential Flow 
About Arbitrary Three-dimensional Bodies,” Douglas Aircraft Co., Ince, 


Report No. E.S. 4062, March 1962. 


Hong, Y.-S. and J.R.e Paulling, "A Procedure far the Computation of Wave- and 


“Motion-induced Forces on Three-dimensional Bodies at Zero Forward Speed,” 


Report to American Bureau of Shipping, May 1979. 


Inglie, ReB. and WeG.e Price, “Calculations of the Yelocity Potential of a 
Translating, Pulsating Source,” Transaction of the Royal Institution of 


Waval Architects, vol. 122, 1980. 


45 


be eet! Davey er waren a 
ed rr 
baw 
r 
a Ai a y ; 
ber nents sonia * i ota “ 
a H 
a my 


10. Abramowitz, M. and L.A. Stegun, “Handbook of Mathematical Functions,” Dover 
Publications, Ince, New York, December 1972. 


ll. Todd, Je, “Evaluation of the Exponential Integral for Large Complex 
Arguments, “Journal of Research of the National Bureau of Standards, vol. 
52. Noe 6, June 1954. 


12. Frank, W., “Oscillation of Cylinders in or below the Free Surface of Deep 
Fluids," David Taylor Naval Ship Research and Development Center Report No. 
2375, October 1967. ; 

13. Brooks, S.He, Ge Weinblum and D.B. Young, “Forces and Moments Experienced 
by a Spheroid Moving Uniformly On or Near The Surface,” David Taylor Naval 
Ship Research and Development Center Report Ns. 384, June 1951. 


14, Nishiyawa, T., “Experimental Investigation of the Effect of Submergence 
Depth upon the Hydrofoil Section Characteristics,” Journal of Society of 
Naval Architects, Japan, vol. 105, 1959. 


15. Wilson, HeB., “Experimental Performance of a Flapped Hydrofoil in Calm 
Water at Low Froude Numbers,” David Taylor Naval Ship Research and 
Devejopment Center Report No. DINSRDC-82/002, June 1982. 


16. Dempsey, E.Me, “Static Stability Characteristics of a Systematic Series of 
Stern Control Surfaces on a Body of Revolution,” David Taylor Naval Ship 
Research and Development Center Report Noe 77-0085, August 1977. 


17. Wadlin, KoL.e and KeWe Christopher, “A Method for Caiculattion of 
Hydrodynamic Lift for Submerged and Planing Rectangular Lifting Surfaces,” 
NACA TR R-14, 1959. 


{8 Newman, J.N-, “The Ferce and Moment on a Slender Body of Revolution Moving 


Near a Wall,” David Taylor Naval Ship Research and Development Center, 
Report No. 2127, December 1965. 


46 


Tre Be Tene ee we me ee em en oe ee eee ewe 


eer! No 
Rr rae bra 


‘bapa aeqna « Ronen jae 0 


oe siaanwel ig hs 
pane ones, A: 


LL bi f 


Soll thy snout secon dunoqethn th tobi 


goto meren ENG eer a they ct bu a tant 


ppl need mp a (aoc ds HPL aR sor Ns Lae Syl lH 
x t } 


"®t 
as 


ATK 
WU el et 


. 


Re 
us 


oe 


Se 
Ras 


xx 
a 
mt 


IND 


SOREN 


ses LS 
eS Ha Pan Nae Pas Co 


Me 
ce 


i Fat ARE Set ENV Nat, 


— 
Ke 


og 


Ue 


TRA 


PIE TEM: 


o's 


as 


fea 


inion 


I sia te i ‘i i 
| Ma Thatiey mse Wen 
M Cll ONG i 


